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Abstract. - We consider a toy model for glassy dynamics of colloidal suspensions: a sin- 
gle Brownian particle diffusing among immobile obstacles. If Gaussian factorization of static 
density fluctuations is assumed, this model can be solved without factorization approximation 
for any dynamic correlation function. The solution differs from that obtained from the ideal 
mode coupling theory (M CT) . The latter is equivalent to including only some, positive definite 
terms in an expression for the memory function. An approximate re-summation of the complete 
expression suggests that, under the assumption of Gaussian factorization of static fluctuations, 
mobile particle's motion is always diffusive. In contrast, MCT predicts that the mobile parti- 
cle becomes localized at a high enough obstacle density. We discuss the implications of these 
results for models for glassy dynamics. 



Introduction. - During the last decade considerable effort has been devoted to simula- 
tional and experimental verification of the mode coupling theory (MCT) of glassy dynamics 
and the glass transition [1-3]. The consensus that emerged from this work is that MCT 
describes in a satisfactory way "weakly" supercooled liquids {i.e. it describes the first few 
decades of slowing down on approaching the glass transition). In particular, MCT has been 
quite successful when applied to concentrated colloidal suspensions [4] , the colloidal glass [5] , 
and gelation [6] transitions. 

Notably, less effort has been devoted to the foundations of the mode coupling theory (see, 
however, Refs. [7-9]). This is somewhat surprising in view of MCT's several well-known 
problems. The most important, fundamental problem is the uncontrolled nature of the basic 
MCT approximation: factorization of a complicated time- dependent pair-density {i.e. four- 
particle) correlation function. 

Recently, we proposed an extension of MCT for dynamics of colloidal suspensions and the 
colloidal glass transition [10]. Our theory includes, in an approximate way, time-dependent 
pair-density fluctuations. It relies upon a factorization approximation that is similar to that 
used in MCT, but is applied at a level of a memory function for the time-dependent pair- 
density correlation function. The theory predicts an ergodicity breaking transition similar to 
that of MCT, but at a higher density. Thus it partially solves another well-known MCT prob- 
lem: overestimation of so-called dynamic feedback effect and the resulting underestimation of 
the colloidal glass transition density. 
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Here, for a simpler, toy model, we go further: we completely avoid using factorization 
approximation for any dynamic correlation function. We only assume Gaussian factorization of 
static correlations [11]. It should be noted that a frequently used approach to glassy dynamics 
is to start from a set of fluctuating hydrodynamics equations which are supplemented by a 
quadratic free energy implying Gaussian static density fluctuations [12,13]. We argue that 
the analysis presented here has implications for such models. 

Since the approach is technically quite involved, we state the main results immediately: 
we derive an essentially exact expression for the time-integrated memory function for a single 
Brownian particle moving among immobile obstacles. We compare this expression with one 
derived from MCT and show that the latter includes a subset of the former's terms: only ex- 
plicitly positive terms from the exact series (i.e. the terms that always increase the effective 
friction felt by the mobile particle) are included within MCT. This is the origin of MCT's 
overestimation of the dynamic feedback effect. An approximate re-summation of the exact 
series suggests that, under the assumption of Gaussian static fluctuations, the mobile parti- 
cle's motion is always diffusive. In contrast, MCT predicts that the mobile particle becomes 
localized at high enough obstacle density. 

This result has important consequences for models used to study glassy dynamics. We 
show here that, if static correlations are Gaussian, a single mobile particle is never localized 
by immobile obstacles. This suggests that a similar fully mobile system (i.e. Gaussian static 
correlations and all particles diffusing) cannot undergo an ergodicity breaking transition. In 
other words, the ergodicity breaking transition predicted for such a system by a mode coupling 
theory is, most probably, an artifact of the factorization approximation. 

Note that this does not mean that MCT is qualitatively wrong for a system with compli- 
cated many-particle static correlations (like, e.g., the hard sphere system). It can be argued 
that terms that cut-off MCT's localization transition (i.e. terms that are neglected in MCT) 
are canceled by other terms that originate from non-Gaussian static correlations. Indeed, 
empirical success of MCT for colloidal systems suggests that this might be the case. It is at 
present unclear how to describe this remarkable cancellation. 

Toy model. We consider one spherical Brownian particle diffusing between N — 1 
immobile, spherically symmetric obstacles. The particle interacts with the obstacles via a 
potential V(r). The obstacles are mechanically identical to the mobile particle. We assume 
that the initial joint probability distribution for the mobile particle and the obstacles is given 
by the equilibrium canonical distribution at temperature T = (fce/3) -1 . The time evolution 
of the system is described by a generalized Smoluchowski equation: 

d d ( d \ 

— P N (Yi\v 2 ,...,Y N ;t) = D — ■ I — -/3Fi J Pjv(ri|r 2 ,...,r iv ;t) 

= OP iv (ri|r 2 ,...,r JV ;t) (1) 

with the initial condition 

PAr(ri|r 2 , ...,r N ;t = 0) = P^(r 1 ,r 2} ...,r N )S(r 1 - r ). (2) 

Here ri denotes the position of the mobile particle and r2,...,r/v denote positions of the 
obstacles. Furthermore, Do is the diffusion coefficient of the mobile particle in the absence of 
the obstacles, and Fi = J2j>i ^lj = ~Sj>i ViV'(ry) i s the force acting on it. Finally, the 
second line in Eq. defines the iV-particle generalized Smoluchowski operator fL 

Our theory starts from the memory function representation of the mobile particle density 
correlation function, Fi(k;t), 

Fi(k;t) = (ni(k)exp(fX)ni(-k)) . (3) 
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Here ni(k) is the Fourier transform of the mobile particle density, rti(k) = e~ Tl , and (. . .) 
denotes the average over the equilibrium probability distribution, P^ 1 . Note the equilibrium 
distribution stands to the right of the quantity being averaged, and all operators act on it as 
well as on everything else. 

To derive the memory function representation we start from an expression for the Laplace 
transform, LT, of the time derivative of F\: 

LT{F\(k) t)) = zF 1 (k,z)-F 1 (k;t = 0) = /m(k)0— !— nx(-k)\ . (4) 



z - n 

Using standard projection operator manipulations we rewrite (0J in the following form: 

LT(F 1 (k;t)) = -D k- (l-iWji(k) J-^ji(-k)Y) • k Fx (*,*)• (5) 

V \ z - QiSiQi / / 

Here ji is a projected current density of the mobile particle, 

j(k) = g 1 ^ (-ik + /3F 1 )e- 4kri , (6) 
Qx = 1 — Pi, and Pi is a projection operator on the mobile particle density subspace, 

A = ■ ■■ nx(-q) )<»i(q) ■■■ = ■■■ ni(-q) )(nx(q) • ■ • . (7) 
q 

Note that in Eq. Q we introduced a summation convention: we sum over all repeated 
wavevectors appearing in adjacent ket, >, and 6ra, <. 

Next we dehne the one-particle irreducible evolution operator, 

n? r = QiD ^Q 1 -(Jt--pF 1 }Q 1 (8) 

and we use the same standard projection operator manipulations to re-write the current- 
current correlation function appearing in Eq. ijSJ in the following form: 

Ji(k, ^w ,( - k, H Ji(k) ^ J,( - k, > (9) 

Combining Eqs. Q and (0 we derive the memory function representation for F\(k; z): 

1 



LT(F 1 (k;t)) = -D k- ^1 + D^ 1 ^ji(k) - _ ji(-k)^ J -kF^k-z) (10) 

Comparing Eq. (|10|) with the standard form of the memory function representation, Fi(k; z) = 
l/(z + 1/(1 + M[ rr (k; z))), we identify the irreducible memory function: 



where k = k/k. 



MT(k; z) = D^k ■ ( j^k)— _j!(-k) ) • k, (11) 
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Next we define n 2 : the part of the joint density of the mobile particle and the obstacles 
that is orthogonal to the mobile particle's density, 

n 2 (qi,q 2 ) =Q 1 J2e~ i ' ll ' ri ~ i ' l2 ' rj - ( 12 ) 
We use the following identity that is exact for two-particle additive interactions: 

ji(-k)) = n 2 (-qi,-q2))(ri2(qi,q2)^2(-q3,-q4)) _1 (n2(q3,q4)ji(-k)), (13) 

where (n. 2 (qi, q 2 )7i 2 (— q3, — q^)^ 1 denotes the kernel of the inverse integral operator. Using 
Eq. l|13|l and an analogous identity for (ji(k) we can express the irreducible memory function 
(1111) in terms of the pair-density correlation function: 

^Mqi,q2;q3,q4;i) = (n 2 (qi,q 2 )exp(fif'Y)n 2 (-q3,-q4)) . (14) 

In turn, this correlation function can be subjected to the same manipulations as the mobile 
particle density correlation function Fi(k;t). Re-tracing the steps between Eqs. © and (fTUf) 
we arrive at the following: 

LT(F 2 (qi,q 2 ;q 3 ,q 4 ;t)) = ~A>qi • ((^(qi, q 2 )n 2 (-q 7 , -qs)) -1 

+ (n 2 (qi,q 2 )n 2 (-q3,-q 4 )} _1 ^j2(q3,q4) z ^^. j^-qs.-qe)^ (15) 

x (n 2 (q5,q6)n 2 (-q7, -qs)) J ■ q7 (n 2 (q7, qs)n 2 (-q9, -qio)> F 2 (q. 9 , qio; q3, q4; z). 
Here j 2 is the two-particle projected current density, 

j 2 (qi,cte) =Q 2 D (JL+pFi) Q^e-**-*-**-*, (16) 
v ri ' j>i 

and Q 2 rr is the two-particle irreducible evolution operator, 

In Eqs. lfTBHi7|) Q 2 = 1 — P% — Pi, with P 1 defined by Eq. JTJl, and P 2 being projection 
operator on n 2 : 

Pi = •■■ra 2 (-qi,-q 2 ) )( n 2(qi,q2)n2(-q3, -q4)} _1 ( n 2 (q 3 ,q4) • • ■ • (18) 



Combining Eqs. (|ll|l . Ijl3(l and l|15|) we can obtain an expression for the memory function 
in terms of the autocorrelation function of the two-particle projected current density: 



Mf (k;z =0)=D o 1 k - (ji(k)n 2 (-qi, -q 2 ))g x 2 qi • (n 2 (qi, q 2 )n 2 (-q 3) -q 4 )) 1 

x ^(n2(q3,q4)^2(-q5,-qe)) + Dq 1 (i 2 (qa , <u) z _ "^ irr J2 (-qs , -qe) \ ^ 

x ("-2(q5,q6)"-2(-q7, -qs)}" 1 • q7<?f 2 (ra 2 (q 7 ,q8)ji(-k)) • k (19) 
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Note that to get expression l|19|) we resorted to a technical approximation that is similar to the 
first Enskog approximation used to approximately invert the Boltzmann collision operator [14] . 
This minor, technical approximation is not required in one dimension where Eq. I|19fl is exact. 

It is clear this procedure can be continued ad infinitum: the two-particle projected current 
correlation function can be expressed in terms of the three-body density correlation function, 
etc. (note that for higher order densities, n m , m > 3, it is advantageous to use ordered 
multiplets of wavevectors, qg < ... < q m )- This is somewhat akin to the well-known continuous 
fraction expansion. It is different from it in that at each step a new function depending on a 
larger number of variables is introduced. 

The resulting expressions simplify greatly if we assume Gaussian factorization of static cor- 
relations. Note that operators Q m remove contributions involving fewer than m + 1 indepen- 
dent connections between two groups of particles [15]. Thus, for example, for the correlations 
of the projected m-particle density we get 

(n m (qi, ...,q m )n m (-k!, -k m )) = N" 1 ^ 1 S{q 2 )...S{q m )S <il , kl ...5 clm ^ m . (20) 

Furthermore, for the density-current correlations of the type (|13|) we get 

(n m (qi, ...,q m )jm-i(-ki, -k m _i)) = inN m ~ 2 (qi - ki) S(q 2 )...S(q m ) 



X! C (*^qi+q^ k l' 5 q2,k2--^q J -l,k J - 1 ^q J + 1 ,k J ---<5q m ,k m _ 1 , (21) 

j=2 



X , 



where n is the number density, n — N/V, and c(q) is the direct correlation function, c(q) = 
(S(q) - l)/(n%)). 

Using the above described procedure, under the assumption of static Gaussian fluctuations, 
we get the following expression for the time- integrated (i.e. z = 0) memory function at k — 

Mr(k = 0; * = 0) = I Y (k . qi) 2 «<«os(«i)*(«o 



qi 



2 \V) 
3! \V) 



2 % , 



1 , f 1 \ 2 c(qi)S(q 1 )c(q 1 )c(q2)S(q2)c(q 2 ) 

qi -^qi • q 2 + k • q 2 ^q 2 • qi 1 



^ 1 qi i2 J |qi + q2p 



^ yk-qi-^qi • q 2 ^— — ( qi +q 2 ) • q 3 + perm 
q1.q2.q3 



X 



1 1 

-qi • 02, — - 

q{ |qi + q2 

c(qi)S(q 1 )c(qi)c(q 2 )S(q2)c(q2)c(q 3 )S(q 3 )c(q 3 ) 
|qi +q2 +qa| 2 



2 



(22) 



where perm, denotes all permutations of the wavevectors' indices. Eq. (122[l is the main 
result of this note. Under the assumption of static Gaussian correlations, this expression is 
essentially exact. The only, technical approximation is the one similar to the first Enskog 
approximation. 

Comparison with MCT. - Expression (1221) should be compared with that resulting from 
the mode coupling theory. The latter approach starts from expression Ijlljl for the memory 
function, rewrites it in terms of the pair density correlation function using (|13|l and JHJ, 
and then resorts to a factorization approximation [16]. This procedure gives the following 
expression for the memory function: 

M^ ICT (k; z) = ^ £ (k • qc(g)) 2 S(q)F 1 (\k -q\;z) (23) 
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Using Eq. (|23[1 we can obtain a series expression for MCT's time-integrated memory function. 
We notice that F^k; z = 0) = (1 + M{ rr {k; z = 0))/(D o k 2 ) and iterate to get 

n S^n x2 c (9i) 5 '(9i) c (ai) 



MS CT (fc = 0;. = 0)4^(k. % 

qi 1 



2 




c{qi)S(qi)c(q 1 )c(q 2 )S(q2)c(q 2 ) 



Iqi + qal 2 



2 



k qi^qi q2, ^ 75 (qi + qa) • qa ) + perm. 

<?r |qi + q2r / / 

x c(qi)S(qi)c{qi)c(q2)S(q2)c(q 2 )c(q 3 )S{q 3 )c(q 3 ) 
|qi+q 2 + q 3 | 2 

where we wrote the resulting series in a form similar to (1221 . It should be noted that from each 
term of 1221) the MCT expression (|24|l includes only these terms that are explicitly positive 
definite. We suggest that the difference between the two expressions is the origin of MCT's 
well-known overestimation of so-called dynamic feedback effect. 

To investigate this last issue further we attempt to re-sum the expression (1221) . To this end 
we resort to an approximation which simplifies angular dependence of individual terms in series 
(|22|l . This approximation will be motivated in detail elsewhere [17]. After the approximation 
we obtain 

2 



V ( k-qi^qi q 2 ...| ■ ^(qi + .. + q m _i) • q m + perm. ] (25) 

c(q 1 )S(qi)c(q 1 )...c{q m )S{q m )c{q m ) 1 

x r— 771 ~ 3^ 2^ c(gi)5(9i)c(gi)...c(g m )5( gm )c(g m ), 

|qi + ... + q m | a qi , q2 ,..., qm 

where d denotes dimensionality of the space. One should note that Eq. (|25|l reproduces 
exactly the first two terms in series ()22|l. Also, Eq. (|25|l is exact in one dimension. 

Using (|25|l we obtain the following simple formula for the mobile particle's diffusion coef- 
ficient: 

D 



D = 



= D exp(-^J2 c ^ S ^)<Qyj (26) 



1 + Mf r (fe = 0;z = 0) 

It is clear that, for a generic interaction potential, the mobile particle's motion is always dif- 
fusive; i.e. the mobile particle is never localized. One should also note a striking resemblance 
between Ij26|l and Deem and Chandler's result [18] for the diffusion coefficient of a single par- 
ticle in Gaussian random media (note that — ksTc(r) is the effective mobile particle-obstacle 
interaction potential [19]). Comparison with numerical simulations of diffusion in Gaussian 
random media showed that the latter result is extremely accurate. 

In contrast, for a generic interaction potential, MCT's equation (|23|l (supplemented by the 
memory function representation for Fi(k; z)) predicts that at sufficiently high obstacle density 
mobile particle's diffusion coefficient vanishes and the particle gets localized. 

It should be emphasized that the absence of the localization transition within the full 
theory and the resulting stark discrepancy with MCT is independent of the approximation 
leading to Eq. H25J1 . The diffusion coefficient in d dimensions is bounded from below by the 
d = 1 diffusion coefficient [20]. In one dimension the series (12211 can be re-summed without 
difficulty. The result is given by formula l)26|l with d = 1. 
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Final remarks. - We showed here that, under the assumption of Gaussian factorization 
of static fluctuations, a single diffusing particle is never mobilized by immobile obstacles, fn 
contrast, MCT predicts localization at a high enough obstacle density. One should notice 
that in some cases, e.g. for a hard sphere diffusing among immobile hard spheres, we expect 
the mobile particle to get localized at sufficiently high obstacle density. Our result suggests 
that this localization is intimately connected to non-Gaussian character of static density fluc- 
tuations. Furthermore, MCT implicitly assumes Gaussian factorization but it neglects terms 
that cut off the localization transition, ft would be of great interest to investigate whether 
in a more accurate theory these terms are canceled by terms originating from non-Gaussian 
character of static density fluctuations. Finally, our results suggest that theories based on 
Gaussian density correlations are not suitable for description of the glass transition. 
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